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Abstract: We study the non-perturbative renormalization group flow of f(R)- 
gravity in three-dimensional Asymptotically Safe Quantum Einstein Gravity. Within 
the conformally reduced approximation, we derive an exact partial differential equa- 
tion governing the RG-scale dependence of the function f(R). This equation is 



shown to possess two isolated and one continuous one-parameter family of scale- 
independent, regular solutions which constitute the natural generalization of RG 
fixed points to the realm of infinite-dimensional theory spaces. All solutions are 
bounded from below and give rise to positive definite kinetic terms. Moreover, they 
admit either one or two UV-relevant deformations, indicating that the correspond- 
ing UV-critical hypersurfaces remain finite dimensional despite the inclusion of an 
infinite number of coupling constants. The impact of our findings on the gravita- 
tional Asymptotic Safety program and its connection to new massive gravity is briefly 
discussed. 
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1. Introduction 

Two central requirements any viable quantum theory of gravity has to meet are 

1. consistency and predictivity at high energies and 

2. the emergence of a classical limit, where the phenomenological successes of 
classical general relativity are recovered. 

Weinberg's Asymptotic Safety scenario [1,2] proposes that these requirements can be 
satisfied within the realm of non-perturbative quantum field theory. In this case, the 
UV-completion of gravity is provided by a non-Gaussian fixed point (NGFP) of the 
renormalization group (RG) flow which renders the theory safe from unphysical diver- 
gences. Provided that the fixed point comes with a finite number of UV-attractive 
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(relevant) directions, this construction is as predictive as a "standard", perturba- 
tively renormalizable quantum field theory. Following the RG-flow emanating from 
the NGFP towards the IR, the second requirement can be fulfilled when the RG- 
trajectories enter into a regime where classical general relativity constitutes a good 
approximation. These ingredients are the foundation of "Quantum Einstein Gravity" 
(QEG), see [3-7] for reviews. 

The main technical tool for investigating Asymptotic Safety is the functional 
renormalization group equation (FRGE) initially constructed for scalar field theory 
[8] and subsequently generalized to gravity [9] 



The FRGE constitutes an exact flow equation governing the change of the effective 
average action (EAA) upon integrating out a small shell of quantum fluctuations 
of the field <3? centered at the IR-cutoff scale k = exp(i). Its solutions determine the 
flow of Tfc on the space of all action functionals. Formally, they interpolate between a 
bare action at the UV-scale k = A and the (standard) effective action (EA) F = Tk=o- 
Notably, applications of ( |1 . 1|) , often known as the Wetterich equation, are not limited 
to gravity, but also play an important role when exploring non-perturbative physics 
in scalar field theory [10-13] and gauge theories [14-16]. 

The key strength of ( |1.1[) is that it allows the construction of non-perturbative 
approximations of the full RG-flow, meaning that they do not involve the expan- 
sion in a small coupling constant. A non-perturbative approximation scheme, which 
is adopted quite frequently, is the truncation of the full to a finite number of 
interaction monomials Oi compatible with the symmetries of the theory 



In the case of gravity, the Oi are often taken as powers of the Ricci scalar R, and 
the truncation order has systematically been increased from the Einstein-Hilbert 
action (N = 1) [9,17-21], to N = 2 [22-24], N = 6 [25,26], and N = 8 [27], up 
to N = 10 [28]. In [29], the operator 3 = R^R^ has been added to the 
N = 2 case and non-local interactions have been considered in [26,30-32]. Moreover, 
the study of quantum effects in the ghost-sector has been initiated in [33], a first 
investigation of "bimetric" flows which discriminate between the background and 
fluctuation fields have been carried out in [34], and RG-flows on surface-terms have 
been considered in [35]. Building on off-diagonal heat-kernel techniques, a general 
strategy for solving (|1 . 1| ) systematically order by order in the derivative expansion 
was proposed in [36, 37] and a physical explanation for Asymptotic Safety based on 
paramagnetic dominance has been put forward in [38]. As a key common result, 
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all these quite sophisticated computations found a suitable NGFP whose associated 
number of relevant directions could be as low as three [7]. Following a similar strategy, 
evidence for a non-trivial fixed point governing the IR-behavior of QEG have been 
reported in [24,39-41]. 

On the technical side, most of these works include a finite number of couplings Oi. 
In this case, one obtains the N + 1 beta-functions (3i, capturing the scale-dependence 
of the Ui, by substituting the ansatz ( |1.2|) into the FRGE and projecting the flow 
on the truncation subspace. The fixed points of this finite-dimensional system are 
given by the zeros of the beta-functions and are thus found by solving the algebraic 
problem (3i(iii) = 0. The UV-relevant deformations of the fixed points are encoded 
in the stability matrix Bij = <9j/3i|FP, obtained by linearizing the RG-flow at a given 
fixed point. In terms of the stability coefficients 9i, defined as minus the eigenvalues 
of B^, the relevant deformations correspond to directions with positive 



The next milestone in studying approximate solutions of ( |1 . 1|) consists in study- 
ing the RG-flow of ^-dependent functions (fRG). In the language of ( |1.2|) this corre- 
sponds to including an infinite number of couplings. At this level, the fixed points 
of the finite-dimensional truncations are promoted to fixed functions which arise as 
k- independent, regular solutions of partial differential equation (PDE) encoding the 
scale-dependence of the functions included in T^. Relevant deformations can then be 
found by studying linear deformations of this PDE around the fixed function. De- 
manding regularity of the deformation leads to a discrete set of relevant operators. 
In scalar field theory this strategy has been first implemented in [8, 42] by deriving 
a non-linear PDE for the scale-dependence of the potential Vk(<ft). In this case the 
PDE includes up to two derivatives of V with respect to <fi. Later on these ideas 
have been refined systematically [43-46] . In [47] the local polynomial approximation 
of the scalar effective average action was used to analyze the dimension-dependent 
existence of universality classes. 

The natural gravitational analogue of the local potential approximation in the 
scalar case is to approximate the gravitational part of the effective average action by 
a generic function of the Ricci scalar 1 



T k [g} = j d d x^f k (R) + . . . . (1.3) 

A crucial difference between the scalar and gravitational theory is that diffeomor- 
phism invariance does not allow to write down a potential for the metric fluctuations, 
i.e., fk{R) doubles as both the kinetic and the potential term of the theory. Substi- 
tuting ( |1.3| ) into ( |1 . 1| ) the flow equation for fk{R) can be obtained by evaluating the 
FRGE on a spherical background, where all possible diffeomorphism invariant inter- 
action monomials can be constructed from the Ricci scalar alone. A PDE governing 



lr The one- loop effective action for /(i?)-gravity has been constructed in [48]. 



- 3 - 



the fc-dependence of fk(R) has been derived in [25,26]. In contrast to the local po- 
tential approximation in the scalar case, this PDE actually contains third derivatives 
of f(R) with respect to R. Using a different cutoff scheme, the asymptotic properties 
of the fixed functions arising in /(i?)-gravity have recently been studied in [49]. 2 The 
complex non-linear nature of the (ordinary) differential equation underlying the fixed 
functions made it very hard to find explicit solutions which are regular everywhere, 
however. 

In this light, this paper studies a somewhat simpler system, the flow of fk(R) in 
conformally reduced gravity in d = 3 spacetime dimensions. In the context of the 
FRGE, conformally reduced flow equations have been studied previously in [50-54]. 
Here it was established that they lead to results quite similar to the full case, where 
the contributions of all metric fluctuations are included. 3 Moreover, working in 
d = 3 considerably simplifies the spectral sum of the Laplacians, so that traces 
of the form Tr[e _sA ] can be evaluated exactly. As it turns out, the "toy model" 
emerging from this setting has all the key features of the PDE governing the scale- 
dependence of fk{R) in full QEG (third order derivatives, non-linearity, fixed and 
moving singularities) while still being simple enough that the combination of analytic 
and numerical methods allows the explicit construction of the corresponding fixed 
functions and their relevant deformations. In this sense, our analysis may be seen 
as an illustrative showcase, providing a comprehensive overview of mechanisms and 
structures expected in the (physically more interesting) case of /(-R)-gravity in four 
dimensions. 

The remainder of our work is organized as follows. The non-linear partial differ- 
ential equation governing the scale-dependence of fk(R) in three-dimensional confor- 
mally reduced QEG is derived in Sect. |[ The scale- independent, regular solutions 
of this flow equation, which constitute the generalization of fixed points to infinite- 
dimensional truncations, are constructed in Sect. |3|. Their properties, including the 
number of UV-relevant deformations, are investigated in Sect. [| and we discuss our 
findings in Sect. ||. 

2. The FRGE on the three-sphere 

We begin the section introducing the basic notions of how to implement a conformally 
reduced gravitational path-integral. The construction will be subsequently adapted 

2 As will be explained in Sect. |^ the key technical difference between [49] and the present work 
is that we do not rely on a summation of the eigenvalues of the Laplacian on a sphere, but use the 
rcsummcd form of heat-kernel [55] . 

3 We expect that the conformal approximation in d = 3 works even better than in d — 4 di- 
mensions, since /(i?)-gravity in the former case contains a single on-shell scalar degree of freedom, 
which is readily captured by this approximation. 
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to the formalism of functional renormalization group. Finally, the flow equation of 
conformally reduced gravity on the three-sphere will be derived. 



2.1 Conformally reduced quantum gravity 

We start our investigation by discussing how conformally reduced gravity can be 
implemented around a definite background geometry. The classical action S griw ['y^ u ] 
for pure gravity is in general a functional of the metric field 7^(2;) . Our aim is 
to construct an effective action where the quantum fluctuations of the metric are 
integrated out. To achieve this in a covariant way, we use the background field 
method that allows to maintain background covariance of the path integral. The 
basic idea of the background field formalism is to split the quantum metric 7 into a 
fixed, but otherwise arbitrary, background field g and a fluctuation field h 

Tfi<w = fjiw + . (2.1) 

It is important to stress that the fluctuation field h^ u is not an infinitesimal pertur- 
bation and hence it needs not to be "small". The fluctuation field will be used as 
the integration variable of the functional integral that defines the quantum theory. 

The theory is said to be "conformally reduced" if we simplify our model by 
restricting the quantum fluctuations to the conformal mode <p = g^h^. In this case 

V - 29i»><P ■ ( 2 - 2 ) 

This is closely related to [56] , where it is suggested that the integration of the confor- 
mal mode should be performed first in a gravitational path integral. An alternative 
procedure for the quantization of the conformal mode involves parametrizing the 
metric with a conformal factor over a fiducial background 

lv» = x 2u g»u (2.3) 

for a certain exponent v. The fluctuationfield x is related to (p via the identification 
X 2u = 1 + \y. It has been shown that it is possible to treat x as the quantum field, 
in fact the quantum theory resulting from this kind of reduction is investigated in 
e.g. [52]. At the quantum level, however, the theories using x an d (p as quantum 
fields are not necessarily equivalent as it was observed in [53]. We will use ip as 
quantum field, since it appears as gauge-invariant scalar mode in the standard York- 
decomposition of the metric. 

The effective action is defined in a well-known way via an integro-differential 
equation 4 



c 



-T[<p;g] _ f ftp e S^ v [gi,u+\9i,uH>\+! x f£(¥>-<^ ( 2 .4) 



4 Here we use the short hand notation f = J d d x y/g. 



- 5 - 



where <fi = (cp) j is the expectation value of (p in presence of the scalar source J = 
We do not need to include ghosts or gauge fixing terms, because diffeomorphism 
invariance is fully broken due to the fixing of the background metric g. In fact, the 
second order variation of any diffeomorphism invariant action, like for example S gr3jV , 
w.r.t. ip generally does not admit any non-trivial zero modes. 

We analyze the RG behavior of the model by constructing a suitable scale de- 
pendent effective average action (EAA). The EAA has a built-in scale dependence 
and can be obtained along the usual lines [9,26]. First we add to the classical action 
in (|2.4j) an infrared (IR) cutoff action, that is quadratic in the dynamical field and 
background-covariant 

AS fc [0; 9] = \f <f>(x)K k [A; g\ <j>{x) . (2.5) 

Here A = — D 2 is the Laplacian operator constructed with the background metric, 
so that background-covariance is always manifest. Additionally, the cutoff operator 
is required to fullfil the boundary-conditions IZk — > for k — > as well as IZk —> k 2 
for k — > oo. This ensures that the EAA interpolates between a "classical" action at 
very high scales and the EA T [0; g] defined in ( |2.4|) at k — 0. 

Infrared and ultraviolet modes are defined w.r.t. the operator A: the role of the 
cutoff term is to effectively suppress the propagation of IR modes, so that UV modes 
are progressively integrated out as the cutoff-scale k decreases. The EAA can then 
be defined through the modified path-integral 



c 



By construction, the EAA has a built-in scale-dependence k and is therefore denoted 



F/C [<j>,g] with the additional label k. If we differentiate the path integral (|2.6|) with 
respect to RG time t = lnfc//c , where k is an arbitrary reference scale, we obtain 
the flow equation 



— Tu\d>: q] — — Tr 
dt L ^' yj 2 



(2.7) 



which is this system's specific form of ( p..l|) . We used the short hand T^' ' [4>;g\ 



to indicate derivatives with respect to the first argument. 
2.2 Computing the variations 

Having the flow equation at our disposal, we now make an ansatz for the EAA. For 
the purpose of the present work, we resort to the single-metric approximation, i.e. 

Y k [<p-g} = Yr\9}, (2-8) 
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with (7) = g = g + hg<f) being the expectation value of the full quantum metric. For 
this EAA we now assume a generic f(R) ansatz 

rf av M = / d d xJgf k {R), (2.9) 

where R is the scalar curvature of the metric g. The usual prefactor (167rGrfc) _1 , 
involving the (scale-dependent) dimensionful Newton's constant G k , is absorbed into 
the now fc-dependent function f k (R). Since the second functional derivative of 1^ 
is needed to compute the flow equation ( |2.7| ), we expand the EAA in a series in h 
around the background field 

T k [g + h] = T k [g,g] + 0{h) + r^tM + 0(h 3 ) . (2.10) 

The conformal reduction is achieved by setting the fluctuations to be h^ u = ^g^cj). 
Performing the second order expansion of Q2.9|) gives 



5 2 T k [g} = J d d x [f k (R) 5 2 ^+f k {R) (2(6Jg)(6R) + ^5 2 R) + y/gfi{R) {5R) 2 ] , 

(2.11) 

where the prime denotes a derivation with respect to the scalar curvature, i.e. 
f' k (R) = The computation is simplified further by choosing a spherically 

symmetric background geometry which e.g. implies D^f^R) = 0. Computing the 
second order expansion in the conformal mode the quadratic part of the series 
Q2.10D becomes 



rr d [<M = I d d xVg~ \<P{x)1C k <P{x) , (2.12) 



with kernel 



** = rl 2 ' 0) [0,^] = ^ [4(d - l) 2 ttD A - 2(d - 1) ((d - 2)f k - Wfi) D 2 

+ {d-2){f k d-ARf k )+AR 2 f h r \ , 

where R is the scalar curvature computed from the background metric g. 

Now it is time to construct the cutoff operator lZ k in such a way that the prop- 
agation of the conformal mode is IR-suppressed. The Hessian (|2.13|) is a function 
of the covariant Laplacian A = —D 2 of the form 

rf 0) = f(-D 2 ; . . . ) , (2.14) 

where all other dependences are denoted by dots for brevity. We choose the cutoff 
implicitly such that the sum of the kernel K and the cutoff operator satisfies 

T {2 ' 0) +TZ k = f(-D 2 + R k ;...). (2.15) 
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Thus lZ k maintains the same tensorial structure as r[ 2 '°' ) , but dresses every Laplacian 
according to 

A-> P k = A + R k . (2.16) 

The function R k = R k (A) is the profile function for the cutoff that contains the 
details of the IR-modes suppression and will be specified later. The function P k = 
Pfc(A) plays the role of a IR-modified propagator for the scalar modes. The cutoff 
action A k S gTllv are determined explicitly by these requirements 



(2.17) 



A fc S grav = ~ J d d x^ Jjj^x) [4(d - Iffi {Pi - D') 
+ 2(d - 1) ((d - 2)f k - 4Rf») R k ] <f>(x) . 
Using fl2.13|) and (|2.17| ) to compute the flow equation (|2.7| ) in the limit = yields 



fr k [0,g] = \Tr[W(A)}, (2.18) 
with Tr [W] = Tr [A/^ 1 ] and 

M =1 [(4(d - lff k ' {P 2 k - D") + 2(d - 1) ((d - 2)/£ - 4i?/£) i? fe )] , (2.19a) 



dt 

.A = [4(d - l) 2 /^P fc 2 + 2(d - 1) ((d - 2)/j[ - 45/0 n 
+ (d-2)(/ fc d-4 J R^)+4 J R 2 /fc] • 



(2.19b) 



It is easy to show that Q2.18 ) coincides with the flow equation for /(i?)-gravity derived 
in [26] in the limit in which the physical scalar is the only dynamical mode of the 
gravitational fluctuations. 

In three spacetime dimensions (d = 3) the function under the trace reads simply 

d 



W(z) = & — , (2.20) 

9kP k + 9kPk + w k K J 

9k = W k , 

~g k = Af k - 16f' k % (2.21) 
w k = AR 2 f k ' -4Rf' k + 3f k . 

2.3 Evaluating the operator traces 

The final step to evaluate the flow equation for the function f k (R) is to perform the 
functional trace of the function W of the operator A appearing the r.h.s. of ( |2.18| ). 
There are in principle different strategies to tackle this task. We briefly outline two 



(g k {Pi -z 2 )+ g k R k 



with coefficients 
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of them that give different results, but only because they ultimately correspond to 
different definitions of how a functional trace should be computed. 

The straightforward possibility is to compute the trace directly by defining it as 
the spectral sum 

Tt[W(A)] = ^W(A B ), (2.22) 

n 

where {A n } are the eigenvalues of the Laplace operator A on the sphere. This has 
been done in [49] applied to the case of four spacetime dimensions (four-sphere). 
However this procedure has its limitations. In particular, if W is a distribution 
(instead of a smooth function), the final result for the trace is a distribution as well 
and therefore it may lead to a flow equation that has some behavior that is hard to 
interpret from a physical point of view and may require further manipulation [49]. 

We decided, instead, to define the trace of any function (and distribution) of the 
Laplacian through the heat-kernel operator e _sA . The trace of a general function of 
A is formally related to the heat-kernel 

oo 

Tr [W(A)] = J ds W{s) Tr [e" sA ] (2.23) 
o 

via the inverse Laplace transform of the function itself defined by 

W{s) = C' 1 [W] (s) . (2.24) 



In contrast to the definition Q2.22 ), it is important to realize that there is no difficulty 



in principle to use the heat-kernel definition of the trace when W is a distribution. In 
fact, any distribution can be seen as the limit of a sequence of functions and therefore 
it is understood that the inverse Laplace transform appearing in ( p.23j ) is the limit 



of the inverse Laplace transforms of the elements of the sequence. It will be evident 
in the following that, in our applications, this limit tends to a well-defined smooth 
function. 

The heat-kernel operator e~ sA occurring in the trace is a well-known object in 
both mathematical and physical literature. Various techniques have been developed 
to compute it, of which we cite only a few [55,57-59]. A standard approach, aiming 
at computing the trace, is to start with an asymptotic expansion for the matrix 
elements of the diagonal heat-kernel 

K{s; x, x) = (x|e" sA |x) = ^ d/2 £ a n (x)s n , (2.25) 



(Atts) 



where the a n are known as deWitt coefficients. If we restrict our attention to a scalar 
theory on a sphere (Euclidean de Sitter space), the asymptotic series (j2.25|) can be 
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resummed up to non-analytic terms. This is achieved by either using a covariant 
expansion of the deWitt coefficients and performing the resummation [55], or by 
means of asymptotic expansions of the Green function at coinciding points [60] from 
which one can read off the expansion of the heat-kernel itself [55]. In d = 3 the exact 
resummation of ( |2.25| ) is particularly simple [55] 

K(s; x, x) = — ^e^ s . (2.26) 



(Airs 



We notice that the result ( p.26|) can also be obtained performing the spectral sum 
( |2.22|) for the case W(z) = e~ sz using the Euler-MacLaurin formula for the summa- 
tion of a series 5 . This observation implies that the summation does not commute 
with the operation of inverse Laplace transform, leading to a crucial difference be- 
tween the definitions (|2~22l) and ( ^23|) . Whether one method is "better" than the 
other is debatable since they are simply different definitions. However, the regularity 
properties inherent in definition ( |2.23| ) have the clear advantage of making the func- 
tional trace fl2.18| ) more regular as a function of the cutoff and therefore less scheme 
dependent. 



Evaluating the trace in (|2.18|) , using the inverse Laplace transform together with 
the heat-kernel (|2.26| ) gives 



Tr [W(A)} 



ds W(s) Tr [e" sA ] 



(47T) 



d 3 x^/g Q 



W 



(2.27) 



where W(x) = W (x — ^) . Moreover Q n denotes the general Mellin-transform of W 
defined by 

1 



Qn [W] 



ds s- n W(s) 



Tin) 



dz z n - x W(z) , 



(2.28) 



where the second equality is true only for n > 0. The proof of ( 2.27|) uses the fact 
that e~ sa represents the translation operator in the Laplace representation of the 
space of functions. In fact, it is easy to prove that 



Q n [W{z + a)] 



ds s- n W{s)e- sa 



(2.29) 



where the notation means that the Mellin-transform is performed with respect to the 
argument z, instead of the full argument of W. 

In our particular case we use the definition (|2.28| ) for the case n — | to evaluate 
the Mellin-transform appearing in ( [2.27|) as 



W 



1 HA R 

— 7"qT" / dz z^W z 

r(f ) h \ 6 



(2.30) 



3 We are grateful to R. Percacci for pointing this out. 
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We stress that up to now the derivation did not rely on any specific property of the 
profile function R k . 

2.4 Cutoff and threshold functions 

Up to now, the profile function R k of the cutoff-operator was kept unspecified. It 
is convenient to choose it to be the so-called optimized cutoff function, which is the 
distribution 



R k ( z ) = (k 2 - z) 6 (k 2 - z) , z = A, (2.31) 

that was originally developed with the purpose of optimizing the convergence of 
the RG-flow [61]. Owed to its technical simplicity, this choice allows many explicit 
exact computations in the context of functional RG and, in particular, of integrals 
like ( |2.30| ). These integrals are often called "threshold functions" since they posses 



a threshold- like structure in the form of a non-trivial denominator. It is useful to 
compute the scale derivative of Rk and the function Pk already at this stage 

^-Rk(z) = 2k 2 6 (k 2 -z), 

dt K J V ; (2.32) 

P k = z + R k = k 2 9 (k 2 - z) + zO (z - k 2 ) . 

Due to the presence of the theta functions, the domain of integration inside ( |2.3C| ) 
reduces to [0, k 2 — c), with c = —R/6. In the particular case of the three-sphere 
R > and therefore the domain [0, k 2 — c] is always non-empty. Using W as defined 
in ( |2.20| ), we obtain 

Q i lw(z+c)]= ml (2 ' 33) 



with coefficients given in Q2.21 ) and 



v k = ~2g k c - g k , (2.34) 
w k = g k k 4 - g k c 2 + 4k 4 g k + g k (k 2 - c) + 2k 2 g k . 

For the optimized cutoff (|2.31| ) the ^-integration can be carried out analytically. We 
denote the basic integrals with a notation similar to Q2.28 ), but take into account 



that the ^-integration is now bounded to the upper limit k 2 — c 

1 r k2 ~ c 1 (k 2 - r) n+m 
Q n [z m ] = / dzz n+m -' = ' l q . (2.35) 



r(n) j r(n) n + m 

Concretely, these integrals are needed in the cases n — | and m — 0, 1, 2 
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2.5 Flow equation 

At this stage we have all the ingredients to write down the full flow equation for the 
/(-R)-truncation. Inserting our ansatz ( |2.9| ) into the l.h.s. of ( |2.18| ), one simply gets 



(2.36) 



where the dot will always denote the derivative w.r.t. the RG time t = log(k/k ). 
Combing our previous results for the functional trace fl2.33p and the definitions ( |2.35| ), 
the r.h.s. of (|2.18|) can be obtained 



1 -Tr[W(A)] = 1 --±- w 

2 2(4^)2 



d 3 x^ Qz [W(x + c)] 



1 



2(4 



7T 2 



d 3 Xv^ 



u k Q3 [z 2 ] + v k Qs [z] + w k Qs [1] 
u k k 4 + v k k 2 + w k 



(2.37) 



Equating ( |2.36 ) and ( ^.37[ ), making use of the explicit form of the integrals (2.35), 
yields the flow equation for f k 



fk 



with 



(6k 2 + Rf 2 ( ai f k + a 2 f> + a 3 f' k + aj f > 
' 567076V (3/ fe + 4 ((k 2 -R)f' k + (-2k 2 + R) 2 f») ) 

ai = 630k 2 , a 2 =2520 (2k 4 - k 2 R) , 

a 3 = 21(6k 2 + R), a 4 =720fc 4 - 432k 2 R - 92R 2 . 



(2.38) 



(2.39) 



This is the desired partial differential equation governing the scale dependence of f k . 

For the purpose of analyzing the flow, we express R and f k (R) in terms of 
dimensionless quantities. Therefore we introduce the dimensionless analogs of R and 

f k (R) 



R = k 2 r, f k (R) = k 3 ip k (R/k 2 ) . 
The derivatives of f k and ip k are related by 

fk = Wk , h = k 3 (ip k + 3<p k - 2r<p' k ) , 



(2.40) 



y k ! 

ft = k'WL 



f k = k (<p' k + </„ - 2r^) , 
ft = W - - 2r^') 



(2.41) 



We can now use ( |2.40| ) and ( 2.41 ) in ( [2.38 ) to write down the complete flow 
equation in terms of dimensionless quantities 



0fc + 3y? fc - 2rip k 



J_ / r\ I cx<p' k + c 2 v k + c 3 (p' k + c 4 (ip'l - 2rtp\ 
v 2 V 6) 3<p k + 4(1 - r)ip' k + 4 (2 - r) 2 ip" k 



(2.42) 
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where the q are polynomials of the dimensionless curvature r 

Cl = j-(36 + r), c 2 =-|-(5r 2 - 234r + 432) , 

f X ° 9 (2.43) 

c 3 = -(6 + r) , C4 =_(6 + r) (30 - 23r) . 

Eq. ( |2.42| ) constitutes the central result of this section. It will be the starting point 
for our search for fixed functions of our functional RG system. 

3. Constructing the fixed functions of the RG-flow 

We now carry out a systematic search for fixed functions of the gravitational RG- 
flow, which arise as regular, scale- independent solutions of eq. ( P-42j ). These solutions 



generalize the fixed points, found on finite dimensional truncations of theory space, 
to the infinite-dimensional theory space of /(-R)-gravity. 

3.1 General structure of the fixed function equation 

Imposing fc-independence of the solution, eq. ( [2.42D reduces to a non-linear ordinary 
differential equation of third order 

y y 7T 2 V 6) 4(2-r)V / + 4(l-r)<£ / + 3^' V ; 

where, in a slight abuse of notation, <p = p(r). The non-linearity of this equation 
makes it very difficult to find its solutions analytically. Thus we resort to a combi- 
nation of analytic and numerical methods to construct its solutions. 

In order that a solution constitutes a valid fixed function, we require that it 
is regular on the entire positive real line r e [0, oof. In terms of the dimensionful 
quantities, this condition reflects the requirement that the fixed function exists for 
all values of k: keeping the dimensionful curvature R of the background fixed, k e 
[0,oo [ i — y r = R/k 2 G ]oo,0]. As we shall see in the following, this regularity 
condition puts severe constraints on the initial data characterizing the solutions of 

We start our search by solving eq. fl3.1| ) for tp'" . Inserting c 4 from eq. ( |2.43| ), we 
obtain 

,„ _ 945A%W,r) 
^ 8(6 + r) (30 - 23r)r ' 1 ' ' 

where the numerator 

M = Uf lin - 7T 2 (1 + r/Q)- 3/2 A/" quad ) (3.3) 
has been decomposed into parts linear and quadratic in ip and its derivatives 

Af quad = (4 (2 - rf ip" + 4{l-r)ip' + 3ip) fa - 2r</j . ^'^ 
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The solutions of (|3.2|) are characterized by three initial conditions for <p, (p f , (p" at 
an initial point r- m i t . This is most easily seen by expressing the solution in a Taylor- 
series expansion around r init . Substituting this expansion into (j3.1|), the higher- 
derivative terms in this expansion are fixed in terms of the three initial conditions 
entering the r.h.s. 

Inspecting the r.h.s. of (|3.2| ), we observe the appearance of "fixed singularities". 
The denominator vanishes at 

30 

r = 0, r = oo, and r = r sing = — « 1.30 . (3.5) 

Expanding the solution at these singular points shows that the poles at r = and r sing 
are of first order while the pole at infinity is second order. These fixed singularities 
motivate to divide the positive real line into two regions 

TZ : < r < r sing , IZ^ : r sing < r < oo , (3.6) 

and use the coordinate x = r _1 on 71^. In addition to these fixed singularities, the 
non-linear nature of ( |3.1|) can give rise to solutions, terminating in "moving singu- 
larities" . In this case, the solution terminates at some r = r term > by developing a 
logarithmic singularity 

ip(r) oc log|r - r tcrm | . (3.7) 

This second type of singular behavior can easily be detected when constructing ex- 
plicit solutions numerically. Based on this discussion of the general structure, we 
now proceed by explicitly constructing the regular solutions of ( |3.1| ) in the next 
subsections. 

3.2 The Gaussian fixed function 

We start by constructing the "trivial" or Gaussian fixed function, which generalizes 
the Gaussian fixed points (GFP) associated with the free or Gaussian theory present 
on any theory space. This "trivial" solution of the fixed function equation is obtained 
as follows. We first introduce a coupling constant c, which controls the interactions 
of the theory. At the level of the function (p(r) this coupling is introduced via the 
rescaling 

f[r(g^)} -> -Mr{g,u)}- (3.8) 

Subsequently, we chose a fixed background and express the metric g^ u through 
the metric fluctuation around this background 

g»u = g» v + \fch iW . (3.9) 

Substituting this expression into the rescaled ip and expanding in h^ u the order 
quadratic in /i„„, encoding the propagator of the theory, is independent of c. Higher 
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order terms, forming the interaction vertices, are proportional to a positive power 
of \fc. Thus, in the limit c — » the action becomes quadratic in h^ Vl so that the 
corresponding path integral is Gaussian. 

When analyzing the limit c — > of the fixed function equation (|3.1|) , we exploit 
that its l.h.s. and r.h.s. are homogeneous of degree one and zero under the rescaling 
Q3.8[ ), respectively. Hence, in the limit of vanishing c the r.h.s., which contains the 
trace-contribution of the Wetterich equation, decouples and one obtains the Gaussian 
fixed point equation 

3(p-2rcp' = 0. (3.10) 

This equation has the solution 

^gfpW =ar 3 / 2 , (3.11) 

where a is a free integration constant. This constant just sets the overall normaliza- 
tion of the free propagator and has no physical significance. 

Owed to the decoupling of the quantum corrections, the result (|3.11|) is easily 
generalized to space-times with dimension d. The numerical coefficients in (|3.10|) just 
encode the canonical mass dimensions of the space-time integral and Ricci-scalar, so 
that the general result can be obtained by substituting 3 — > d in the first term. Thus 
the (i-dimensional Gaussian fixed point solution is given by 

l p GFP (r) = ar d / 2 . (3.12) 

For d = 4, this general result reduces to the Gaussian fixed point solution found in 
[49]. Notably, the solutions ( |3.12| ) are distinguished by the property that the dimen- 
sionful coupling constants multiplying the R d / 2 -term are powercounting marginal, 
i.e., of mass- dimension zero. This is expected from the fact that, for the free theory 
at the Gaussian fixed point, there are no quantum corrections that could lead to an 
anomalous dimension of the operators. 

3.3 Non-Gaussian fixed functions 

The remainder of this section is devoted to the search of non-trivial regular solutions 
of eq. ( |3.1| ). The analysis of the singularity structure ( |3.5| ) showed that the fixed 
function equation possesses a double pole at r = oo, while the pole at r = is of 
first order. We thus expect that imposing the regularity condition at large r will fix 
two of the three initial conditions characterizing the solutions, while an analogous 
analysis at small r may leave two parameters undetermined. Hence we will start in 
the IR (r ^> 1) and continue our analysis towards the UV (r <C 1). 

Asymptotic solutions at large r 

We start with fixing the asymptotic behavior of the fixed functions for r — > oo. For 
this purpose, we make the ansatz 

(p~cr a (l + u 1 r- 1 + ...), (3.13) 
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where the exponent a determines the asymptotic behavior of the solution, and the 
dots denote subleading terms. Substituting this ansatz into (|3.1| ), the asymptotic 
expansion yields 

cr a (3-2a) = I 5670 ^ 2 M (3 . 14) 

945-67«l 3/2 , _ o In 

y 34020 VS* 3 r + • • • a — a/z. 

When expanding the r.h.s., the limit a = 3/2 is non-trivial, since in this case the 
leading terms in the numerator and denominator both vanish, such that the coeffi- 
cient U\ appearing at subleading order actually enters into the equation. Comparing 
powers of r, ( J3.14 ) establishes that for a > 3/2 the only solution of this equation is 



c = 0. Thus the highest power of r allowed in the asymptotic expansion is a = 3/2. 
This case allows the construction of a non-trivial solution (c ^ 0) by adjusting the 
subleading coefficient u\. Notably, the asymptotic behavior r 3 / 2 reproduces the large 
R limit of the one-loop effective action in three dimensions [55]. 

Together with the square-root appearing in the fixed function equation, this 
asymptotic behavior motivates redefining 

cp(r) = (1 +r/6) 3/2 y(r) . (3.15) 



As it turns out, the prefactor in ( |3.15| ) captures the asymptotic behavior of ip( 



r 



for both small and large values of r. The profile function y(r) then interpolates 
continuously between these two asymptotic regimes. 

In order to proceed with the analysis, we adapt the differential equation for the 
profile function y(r) to the patch 7^oo by changing variables x = 1/r. The function 
v(x) = y(l/x) then satisfies the non-linear third order equation 

'" = N{v,v',v",x) 
V 240 x 2 (1 + 6x) 4 (x - x sing ) ' 1 ' ' 

Here the prime denotes a derivative with respect to the argument x, x sing = l/r sing = 
23/30, and the numerator is given by 



M = c (18 (19 - 558a;) v - c c l v' + 2 x c 2 c 2 v") 

- 3780tt 2 x (9v + c v') ((3 - 78x) v - 2x 5 c 3 v - 2x c 2 (1 - 2x) 2 v") . 

The coefficients Sj are polynomials in x 

d = 1 + 6x , ci = 67 - 1434 x - 72360 x 2 + 103680 x 3 , 

c 2 = 113 + 4302 x - 6480x 2 , c 3 = 21 - 58x + 48x 2 . 



(3.17) 



(3.18) 



As already anticipated, the denominator of ( |3.16| ) has a double zero at x 



0. Insisting that the asymptotic behavior of the solution is given by ( |3.15 ), this 
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divergence fixes two of the three initial conditions characterizing the solutions of the 
fixed function equation. For this purpose we expand v(x) around x = 

v (x) = c(l + v 1 x + v 2 x 2 ) +0{x 3 ). (3.19) 



Substituting this ansatz into fl3.16| ) and expanding the resulting equation at x = 
results in a Laurent series starting with x~ 2 . Requiring that the series coefficients of 
x~ 2 and x^ 1 vanish fixes the unknown coefficients Vi, v% in terms of c 

342 90 

vi = -^r , v 2 = — (376 + 19845 vr 2 c) . (3.20) 

The parameter c remains undetermined at this stage. In the next subsection, this 
parameter will be constrained by demanding regularity of the solutions at r sing and 
r = 0. 

Regularity at r = r s i ng 

At this stage, we have restricted the set of candidate fixed functions to a one- 
parameter family of functions v(x;c), depending on the initial condition c imple- 
mented at x = 0. In the next step, we apply a numerical shooting method to extend 
the asymptotic solutions (|3.19|) to the interval x e [0,x S i ng ]. 



Owed to the first order pole in fl3.16|) , we a priori expect that the solution for a 
generic value c will become singular at x — x sing . Expanding the solution in a Taylor 

v(x; c) = — a n {c) (x - x sing ) w , (3.21) 

n=0 U ' 

and substituting this ansatz into ( |3.16| ), the l.h.s. of the differential equation remains 
regular while the r.h.s. generically has a first order pole. The cancellation of the 
residue then gives rise to an additional constraint on the admissible fixed functions. 
The shooting method allows to phrase this constraint in terms of the free parameter 
c, fixed at x — 0. Based on these preliminary considerations, we expect that only a 
discrete set of the candidate solutions is able to pass a; S i ng without becoming singular, 
thus limiting the allowed values of c to a discrete set. 

In practice, we construct v(x;c) by numerically integrating Q3.16] ) up to x ter m = 
^smg — applying a BDF-algorithm. Varying c, this gives rise to a one-parameter 
family of solutions v(x;c), defined on the interval x G [0,x term ]. Substituting these 
solutions into ( 3-171 ), we define 

6(c) = lim J\f(v(x; c), v'(x; c), v"(x; c); x) . (3.22) 

X ^3?sing 

The vanishing of the residue at x — x s - mg is implied by the zeros of 6(c). In other 
words, solutions v(x;c), which are regular at x S i ng correspond to critical values c cr i t 
where 5(c crit ) = 0. 
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Figure 1: The function 6(c), eq. (3.22). Zeros of 5(c) correspond to solutions v(x;c) 



which pass through the singularity at r 



' sing 



and can be continued to the patch TZq. For 



c > there are three isolated zeros located at ci )Cr i t = 3.90 x 10 , C2, C rit = 4.41 x 10 
and C3 )Cr it = 0.255. In addition there is one continuous window when 0.648 < c < 2.668. 



The function 5(c) constructed via this algorithm is shown in Fig. [I]. Here the 
sampling of 5(c) has been carried out in steps of Ac = 10~ 4 . The figure displays 
three isolated zeros at 

c ljCrit = 3.90 x If!" 4 , c 2 , crit = 4.41 x 10" 2 , c 3 , CTi t = 0.255 , (3.23) 

and a continuous window of zeros located in the interval 

c-, c rit <c< c +iCrit , c_ crit = 0.648 , c +jCrit = 2.668 . (3.24) 

We also investigated the dependence of the critical values Cj jCrit on the parameter 
e and found that their position is stable, if e is chosen sufficiently small. When 
decreasing e, the function 5(c) is becoming steeper at the isolated critical points, 
while inside the window it is continuously approaching zero. 

While the existence of isolated solutions has already been anticipated in the 
discussion above, the appearance of the continuous window (|3.24j) is somewhat sur- 



prising. In this case, 5(c) = acts as an attractor for an entire interval of initial 
conditions. Since M is build from v (x) and its first and second derivatives, which 
balance at x s ; ng the solutions v (x; c) in this interval are physically different also when 
continued to TZq, i.e., one does not "lose memory" by passing through this attractor. 

At this stage, the following remarks are in order. We have also investigated the 
behavior of 5(c) for negative c. In this regime ( [3.22 ) is positive definite, implying 



that these solutions cannot be extended into the patch TZq. Moreover, by expanding 
( p.42|) for small values c, which essentially corresponds to linearizing ( |3.16|) , the 



solutions v(x; c) can be found analytically. Substituting the analytic expression into 
5(c) one can show that for c < ci jCr i t there is only the trivial zero c = 0. In this 
regime the analytic and numerical studies match continuously, ensuring that there 
is no accumulation of zeros for c C 1. 
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0.26 ■ 
0.25 ■ 
0.24 ■ 
0.23 ■ 
0.22 ■ 
0.21 ■ 

1.26 1.28 1.30 1.32 1.34 1.36 1.38 1.40 1 

Figure 2: Numerical continuation of the candidate fixed, function C3 ;Cr it = 0.255 to the 
IR-patch TZq . The solution terminates in a moving singularity at rt er m = 1-26 and cannot 
be continued to r = 0. 




Continuation to r = and moving singularities 



In the final step, the candidate fixed functions encoded in ( 3.23Q and ( 3.24j) are 



continued to the UV-patch IZo- There are two potential obstacles, which may prevent 
the solutions to constitute valid fixed functions: 

1. moving singularities at r tcrm > 0. 

2. the first order pole of ( |3.2| ) at r = 0. 

Following the strategy of the previous subsection, we investigate these final condi- 



tions by numerically integrating ( |3.2|) with initial conditions at r sing , employing an 
algorithm that takes the stiffness of the ODE at r = into account. Practically, 
we perform the continuation of v(x\c cr - lt ) to TZq by first determining w(c;c cr i t ) and 
its first and second derivative at x sing — e. These values are then used as initial 
conditions for the profile function y(r) at r sing — e, taking the proper Jacobians into 
account. Typically, we work with e = 10 -7 , and we have checked that our results are 
insensitive to e as long as it is sufficiently (but not too) small. 

The only case in which we encounter a moving singularity is the solution ema- 
nating from C3 iCr i t . The corresponding numerical integration is shown in Fig. ^|. This 
solution develops a logarithmic singularity of the form ( |3.7| ) at r term = 1.26 and can 
therefore not be completed in a regular solution extending down to r = 0. Notably, 
the other candidate fixed functions ( |3.23| ) and in particular all solutions from the 
continuous window ( |3.24| ) are free from this type of singularities. 

Finally, we encounter the fixed singularity at r = 0. In analogy to the singu- 
larity at r sing , we require that this pole in ( |3.1|) is compensated by a zero in the 
numerator (|3.3|) . We then substitute our numerical solutions into J\f(tp,(p',tp",r). 
For the isolated solutions ci jCr i t and C2, C rit and the two boundaries of the continuous 
window, the result is shown in the left and the right panel of Fig. [| respectively. At 
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0.2 0.4 OS CTS- 



Figure 3: The function J\f evaluated along the critical solutions c\ cr jt (black curve) and 
C2,crit (red curve) and the boundaries of the continuous window c_ iC rit (black curve) and 
c+,crit (red curve) are shown in the left and right panel, respectively. For presentational 
purposes, M(r; ci. cr it) has been magnified by a factor 10 3 . The vanishing of the numerator 
at r = ensures that the solutions are regular on the entire interval ^ v <C ^sine* 

r = 0, J\f(<p, <p', <p", r; Cj jCr ;t) vanishes, establishing that the last condition is satisfied 
automatically for all remaining fixed functions. 



In summary, we have established that besides the Gaussian fixed function ( |3.11 



the non-linear fixed function equation ( |3.1| ) gives rise to two isolated and one contin- 
uous family of fixed functions which are regular on the entire positive real axis r > 0. 
They can conveniently be characterized by their asymptotics in the IR (r ^> 1): 



Cl crit = 3.90 x 10" 



C 2 ,crit = 4.41 X 10" 



(3.25) 



and 



c- CT it < c < c +iCrit , c_ crit = 0.648 , c +jCrit = 2.668 . (3.26) 

The existence of these fixed functions is the central result of this section and we will 
continue by studying their properties in the remainder of this work. 



4. Properties of the fixed functions 



As the central result obtained so far, we established that the flow equation (|2.42j ) of 
conformally reduced QEG in d = 3 possesses two isolated and one continuum family 
of regular fixed functions. The characteristic properties of these solutions will be 
discussed in this section. 



4.1 Characteristics of the numerical solutions 

We start by displaying the regular fixed functions <fi,*(r) arising from the critical 
values c i cri t, eqs. ( |3.25| ) and ( |3.26|) , in Fig. |[ The corresponding profile functions y(r) 
introduced in ( |3.15| ) are shown in the right diagrams. They interpolate continuously 
between the UV (r 1) and the IR (r ^> l)-regime 6 , where they become constant. 

6 In a slight abuse of language we will refer to the expansions of our fixed functions for large and 
small values k as UV and IR. The corresponding fixed functional T* is of course independent of k. 
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Figure 4: All complete, regular fixed functions ( 3.25| ) and (|3.26| ). The upper panels 
show the isolated fixed functions characterized by ci. cr it (black) and C2 iCr it (red). The lower 
panels display the boundary solutions of the continuous window c_ jCr ;t (black) and c+ jCr ;t 
(red). For presentational purposes the ci jCr ;t-solution has been multiplied by a factor 100. 



The crossover occurs when the RG-scale is of the same order of magnitude as the 
background curvature, i.e., r = R/k 2 ps 1. 

The fixed functions ipi t * (r) resulting from these profile functions are shown in the 
left diagrams of Fig. @. They are all bounded from below and increase monotonically 
with increasing r. The asymptotic behavior of the profile functions furthermore 
implies that the <fi,*(r) interpolate continuously between 

^(r)~u* , r«l (UV) 

^(r)~6- 3 / 2 c, crit r 3 / 2 , r»l (IR) . 



Here the IR-coefficients Cj, C rit are given in eqs. ( |3.25 ) and ( |3.26| ), while the constants 
u* Q determining the UV-asymptotics are listed in the first line of Table [TJ. When com- 
paring the four solutions, there is an obvious qualitative difference between (pi t *(r) 
and the other fixed functions: <£>i,*(r) is singled out by the fact that it is the only 
solution that remains positive in the UV. Thus it is the only fixed function that 
remains positive definite throughout. 

At this stage it is illustrative to restore the dimensionful quantities and investi- 
gate the asymptotic behavior of the conformal field theories associated with the fixed 
functions in the IR. Substituting the IR-asymptotics ( |4.1| ) into ( j2.4U| ) and restoring 
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the spacetime integral via ( |2.9|) , the effective action T becomes 

r * = ^2 / d3x V9Ci,cntR 3/2 ■ (4.2) 

The effective action resulting from our fixed functions is precisely given by the power- 
counting marginal operator of the theory with the coupling constant determined by 
Cj, C rit- Notably, it has the same form as the one-loop effective action [55]. A similar 
structure for the IR-asymptotics of r* has also been observed in four dimensions [49] . 
In this case it was, however, not possible to obtain the allowed values Cj jCr it since the 
complete regular solutions of the differential equation describing the fixed points has 
not been constructed. 

We stress that none of the properties discussed here are put in "by hand" . They 
all arise from solving the fixed function equation ( |3.1|) and this have the status of 
predictions from the quantum theory. 

4.2 Asymptotic expansions 

Refs. [25-27,62] studied the RG-ffow of -gravity by implementing a polynomial 
ansatz for the function f(R) and subsequently analyzing the properties of the result- 
ing NGFP. 7 It is thus instructive to compare the full numerical solutions shown in 
Fig. [| and their polynomial approximation found by solving the fixed point equation 
recursively. 

Polynomial expansion for small r 

In order to capture the asymptotic behavior of the fixed functions in the UV, we 
expand the fixed function in a Taylor series at r = up to a fixed order N 

N 
n=0 

Substituting this ansatz, the expansion of (|3.1| ) at r = gives rise to a first order pole. 
Following the discussion of Section |3.3|, this pole fixes one of the three free parameters 



characterizing the fixed functions. The remaining two constants of integration can 
then be taken as u* and u\. By solving the expanded fixed function equation order 
by order in r, the coefficients u n for n > 2 can then be found recursively as functions 
of these free parameters, i.e., the coefficient u* 2 is fixed by a condition at order r°, etc. 
It thereby turns out that the resulting system of equations has a triangular structure: 
the first equation is linear in u\ and thus gives a unique value for u\ in terms of Uq 
and it*. The next order equation is again linear in u^, and uniquely determines 
in terms of the lower order couplings, and so on. This structure guarantees that the 
solutions are characterized uniquely by Uq and u*. 



7 For a similar study in scalar-tensor theories see [63]. 
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c l,crit 


C2,crit 


C— ,crit 


C+,crit 


u* 


1.185 X 1(T 3 


-2.635 x 10~ 2 


-5.099 x 10- 1 


-2.056 


u\ 


8.466 x 10" 5 


6.917 x 10" 2 


8.211 x lO" 1 


2.996 


^2 fit 
U 2,rec 


8.593 x 10" 6 
8.595 x 10" 6 


-5.697 x 10" 3 
-6.143 x lO" 3 


1.349 
-5.515 x 10~ 2 


10.554 
-0.182 


^3, fit 
n 3,rec 


-1.080 x 10~ 6 
-1.098 x 10~ 6 


-3.922 x lO" 3 
-8.828 x 10~ 4 


-8.109 
-5.389 x 10- 3 


-56.701 
-1.450 x 10~ 2 


U 4,fit 
M 4,rec 


2.944 x 10~ 7 
4.104 x 10~ 7 


9.770 x 10- 3 
3.246 x 10~ 4 


22.874 
-1.549 x 10~ 3 


149.672 
-3.549 x 10~ 3 


U 5,fit 
11* 


-1.222 x 10~ 8 
2.474 x 10~ 8 


-1.874 x 10~ 2 
-1.566 x 10~ 4 


-37.364 
-6.187 x 10~ 4 


-230.235 
-1.2289 x 10" 3 



Table 1: Series coefficients appearing in the expansion ( |4.3| ) of the fixed functions ipi^(r). 
The two free parameters u* and u*^ characterize the fixed functions in the UV-regime 
for r<l. Together with the higher order momenta u* fit they are obtained by fitting the 
polynomial (4.3) to the numerical solutions 99^* (r), while the values u* nrec are determined 
recursively by solving the system of linear equations arising from expanding ( |3.1|) up to 
order r n_2 . 



With the complete numerical solutions at our disposal, we can fit the polyno- 
mials ( |4.3| ) (with N = 50) to the fixed functions (r) restricted to the UV-patch 
TZq . The parameters u* Q and u±, which characterize the solutions, and the first four 
subsequent coefficients u* m are shown in Table [I]. Taking the parameters u* Q and 
u* x (the first index labels the various solutions Cj )Cr it) as initial conditions for the 
polynomially expanded fixed function equation, we also determine the higher order 
expansion coefficients recursively. These are summerized by the entries u* rec of Ta- 
ble |l| Comparing the expansions obtained in these two ways, we observe that the 
expansion coefficients of <£>i,*(r) show a reasonable agreement, while the expansions 
of the other solutions already show some level of disagreement at order u* 2 . 

These observations can be formalized by studying the convergence radius of the 
recursive solutions r con . We determine r con via the quotient criterion 

r con = lim . (4.4) 

For the fixed function </?i,*(r) the expansion coefficients u n are easily obtained up to 
order 200. The resulting first 200 terms in the series (£OJ) are then shown in the left 
diagram of Fig. || For n > 50 the series converges rapidly. Fitting a Laurent series 
including a constant, r _1 and r~ 2 term, we obtain 

r con = 1.30707 ^r sing . (4.5) 

Thus the polynomial expansion for ip\^(r) converges up to the singularity r S i ng . This 
result is also demonstrated in the right panel of Fig. [5], which displays the direct com- 



- 23 - 



0.001 40 r 




50 100 150 20(1 



Figure 5: Series coefficients ( |4.4[) , arising from the recursive solution of the fixed function 
equation with initial conditions u* and u\ 1 . The series converges up to r con = r s i ng . This 
is confirmed in the right panel which compares (p± *(r) obtained numerically (black) and 
recursively (red). 




Figure 6: Comparison between the fixed function f2,*( r ) obtained numerically (black) 
and its polynomial approximation ( |4.3| ) truncated at order N = 5 (red), N = 10 (blue), 
and N = 15 (magenta). Notably the quality of the polynomial approximation decreases 
with increasing order of the approximation, indicating that its radius of convergence j| ) 
is actually zero. 



parison between the fixed function obtained numerically (black) and the polynomial 
approximation (red). 

Applying the quotient criterion to the other fixed functions quickly reveals that 
their radius of convergence is actually zero. We exemplify this result in Fig. 
which compares the numerical solution v 9 2,*( r ) and the corresponding polynomial 
approximation with initial values u^q and u^i- Increasing the order of the Taylor- 
expansion systematically actually decreases the quality of the approximation. This 
behavior is prototypical for the non-convergent nature of the series. 



Physical coupling constants in the UV 

The UV-behavior of the fixed functions is governed by the expansion ( |4.3| ) around 
r = 0. Restoring dimensionful quantities, the UV-asymptotics of the average action 
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Q,crit 


A* 


9* 


9 

mi 


C-ljCrit 


3.90 x 10" 4 


-7.00 


-235 


1.23 


C2,crit 


4.41 x 10" 2 


0.19 


-0.29 


-8.6 x lO" 3 


C— ,crit 


0.648 


0.31 


-2.4 x 10" 2 


7.6 x 10~ 2 


C+,crit 


2.668 


0.34 


-6.6 x 10" 3 


3.5 x 10^ 2 



Table 2: Physical couplings arising from the polynomial expansion of the fixed functions 
<Pi * (r) in the UV. The values are obtained by substituting the coefficients it* fit into the 
relations ([0| and ( p~lc| ). 



r*[<?] resulting from the fixed functions fi t *(r) has the curvature expansion 

1 



lim T*[g] 

fc— >-oo 



16ttG 



-R + 2A k ) + u 2 R A + ... 



N 



(4.6) 



where A, Gn and u 2 are the dimensionful cosmological constant, Newtons constant 
and _R 2 -coupling, respectively. They are related to their dimensionless counterparts 
by 

A fc = X k k 2 , G k = g k k~ l , u 2 = u 2 k . (4.7) 

A comparison between the polynomial expansion ( |4.3| ) and the Einstein-Hilbert ac- 
tion allows to express the dimensionless physical couplings to the expansion coeffi- 
cients 



A* 



Ur 



1 



9* = -TTT-r • (4-8) 

2u\ 1&TTU\ 

Notably, the structure of ( |4.6|) closely resembles the new massive gravity [64, 65]. 8 
In order to fully exploit this connection, it is useful to trade the coupling u 2 for the 
mass of the massive gravity mode. This relation is found as follows. Substituting 
the f k (R) resulting from the expansion ( f4.6| ) into (|2.13 ) and afterwards taking the 
flat-space limit by setting R = 0, A = 0, we obtain 



/CI 



flat space 



2(d 



d? 



— 2/2 

-u 2 p [p 



(d-2) ui 
4(d-l) ti~2 



(4.9) 



Performing a Wick-rotation the zero in (|4.9|) yields the mass of the "massive graviton" 
in d = 3 



VTLr 



■k 2 = m 2 k 2 



(4.10) 



8u 2 Su 2 

The values of the physical dimensionless couplings ( |4.8|) and (|4.10|) arising from 
the polynomial expansion of our fixed functions in the UV are collected in the third, 
fourth, and fifth column of Table [| In addition their values across the continuous 
window of fixed functions ( ft -26] ) are displayed in Fig. [7|. 



3 For a treatment in the framework of the FRGE see [66]. 
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Figure 7: Value of the dimensionless cosmological constant A* (left), Newtons constant 
(center), and the massive graviton mass ml (right), arising from the UV-expansion of 
the fixed functions <fi*(r) in the continuous window ( |3.26| ). 



We find that the the mass of the massive graviton is always positive while New- 
tons constant remains negative throughout. While the later result looks unphysical 
at first sight, one should keep in mind though that the kinetic term of the conformal 
factor (the latter driving all the quantum effects considered in this paper) actually 
comes with the "wrong sign" for a kinetic term. This is precisely compensated by 
the negative value of Newtons constant, so that the fixed function produces a kinetic 
term that is positive definite. 

4.3 Linear perturbations and relevant operators 

An additional characteristic of the fixed functions ip* (r) are their UV-relevant defor- 
mations, which will be studied in this subsection. For this purpose, we return to the 
full, t-dependent partial differential equation (|2.42p and substitute the ansatz 



<p(t,r) = <p*(r) + ee~ yt v(r) . (4.11) 

Here 9 denotes the stability coefficient associated with the perturbation v (r), y?*(r) 
specifies the fixed function for which we study the deformations, and e is a book- 
keeping parameter, indicating that the deformations under consideration are actual 
infinitesimal. By definition, UV-relevant deformations approach y2*(r) as t — > oo and 
are thus characterized by 9 > 0. 

Substituting the ansatz ( [4.1 1|) into (|2.42|) and expanding in e, the fixed function 



equation ( |3.1|) is recovered for e = 0. The information about relevant deformations 
is encoded in the linear third order differential equation for v(r) obtained at order e: 



(3 - 9)v - 2rv> - 1 (1 + M 3/2 CW + CW ' - 9cW - C4 {9V " + ^ 
[6 9)v Zrv - - (1 + g j 3^ + 4( i_ r )^ + 4 ( 2_ r )2^ 

(3</?* - 2?VJ (3v + 4(1 - r)v' + 4(2 - r) V) 
3^ + 4(l-r) ^ + 4(2-r)V: ' 



(4.12) 



Here the coefficients q are given in (|2.43|) and we used that <p* satisfies the fixed 
function equation ( |3.1| ) in order to simplify the resulting formula. Eq. ( |4.12| ) depends 
explicitly on the fixed function tp* and parametrically on the stability coefficient 9. 
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Cl,crit 


c 2,crit 


C— ,crit 


C+,crit 




4.896 








02 


0.287 


0.058 


0.006 


0.002 



Table 3: The positive critical exponents associated with the UV-relevant deformations 
of the fixed functions shown in Fig. ||. 



Analogously to the study of fixed functions, relevant deformations v(r) are given by 
solutions of ( |4.12|) which are regular on the entire interval r £ [0, oof. We expect 
that this regularity condition puts a constraint on the admissible values 0, so that 
there is only a finite number of relevant deformations. 

When searching for the values which give rise to regular solutions we proceed in 



complete analogy to Section 3.3. Following the strategy for the scalar case outlined 



in [42], we first study the asymptotic behavior of ( [4.12 ) for large r. In this limit 



we identify the three fundamental solutions of the linear differential equation, which 
scale as 

v(r) ~ r° , v{r) ~ r 113 / 92 , and v(r) ~ r^' 2 , (4.13) 

respectively. The last solution is the asymptotic behavior expected from scaling 
arguments [42]. We thus fix the asymptotic behavior of the deformations according 
to v(r) ~ r^ 3 ~ 9 ^ 2 and set the coefficients for the other asymptotic solutions to zero. 
After fixing the asymptotics of our solutions we integrate ( 4.12| ) numerically and 



single out the parameters for which the solution remains regular at x — x S i ng . 
Restricting ourselves to the case of relevant deformations with 9 > the result is 
summarized in Table |[ We find that the fixed function (pi^(r) admits two relevant 
deformations while for the other fixed functions there is only one admissible value 
9 > 0. Thus the regularity condition indeed restricts the allowed values of to a 
discrete and finite number. This central result constitutes strong evidence, that the 
UV-critical surface of the corresponding fixed functionals is actually finite dimen- 
sional, despite the inclusion of an infinite dimensional set of coupling constants in 
our truncation ansatz. 

At this stage, the following remark is in order. The ansatz (£ d]) is restricted to 



real stability coefficients. In order to encode complex stability coefficients, eq. ( [4.11| ) 
needs to be generalized to complex deformations v(r) = fR C (r) + ivi m (r) according 
to 

(p(t, r) = <p,(r) + \e (e~ {e ' +te "^ v(r) + c.c.) . (4.14) 



Substituting this generalized ansatz into Q2.42 ) results in a coupled system of third 



order differential equations which depend on the real and imaginary part of v(r) and 
two parameters 0' and 0". We investigated this system numerically, but were not 
able to find any regular deformations with complex critical exponents 9' + %9" . 
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5. Conclusions and outlook 



In this work we used the exact functional RG equation [9] to derive a non-linear par- 
tial differential equation (PDE) governing the RG-flow of -gravity in conformally 
reduced Quantum Einstein Gravity (QEG) in d — 3 dimensions. In comparison to 
the previous works [26,49], our main technical improvement is the use of the exact 
heat-kernel formula on the sphere [55] which allows the construction of this PDE 
without exploiting specific properties of the regulator function IZk- 

Our key interest is in scale-independent, regular solutions of this PDE which 
in this particular case simplifies to an ordinary non-linear third order differential 
equation. The resulting fixed functions constitute the generalization of the fixed 
points, known from finite-dimensional truncations of the effective average action, to 
the realm of infinitely many coupling constants. By combining analytic and numeric 
methods, we construct two isolated and one continuous one-parameter family of 
fixed functions. These constitute the first non-trivial fixed functions found in the 
gravitational framework. The effective action associated with the fixed functions has 
the form 

r*[#] = Q- 3/2 J d 3 x^c hCvit R 3/2 } (5.1) 



with the critical values Cj )Cr i t given in ( |3.25| ) and ( p.26| ) . Notably, this form of the 



effective action is the same as for the one-loop effective action found previously [55]. 

For values of the RG-scale k much larger than the curvature R, the solutions 
admit a polynomial expansion of the form f(R) = Xl^Lo ■ Based on this ex- 
pansion, all solutions give rise to a negative Newtons constant, which is required for 
the positivity of the conformal factor kinetic term. Moreover, one of the discrete 
and all continuous fixed point solution have a positive massive graviton mass. By 
studying linear deformations around the fixed function, we furthermore establish that 
they come with a finite number of UV-relevant deformations. Table |3| indicates that 
there are two and one relevant deformations of these solutions respectively. Thus 
the UV-critical surface of the fixed functions is actually finite- dimensional, despite 
the inclusion of an infinite number of coupling constants in the truncation. This 
result further strengthens the evidence that Asympotically Safe Quantum Gravity 
is actually a predictive quantum theory of gravity with a small number of relevant 
coupling constants. 

The physically most interesting properties are shown by the fixed function char- 
acterized by Ci jCri t- In addition to the properties found above, the corresponding 
fixed function is positive definite with a global minimum at zero curvature. More- 
over, it can be shown that the polynomial expansion for RG-scales much larger 
than the curvature actually has a finite radius of convergence determined by the 
singularity structure of the underlying differential equation. Thus the polynomial 
/(/^-approximation, pioneered in [25] and subsequently refined in [26-28], may give 
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A* 


9* 


2 

mi 


6i 9 2 


^3 


NGFP 
UNGFPi 
UNGFP 2 


0.019 
0.364 
0.099 


0.188 
0.147 
0.216 


-1.67 
-0.31 
-0.22 


8.39 1.86 

1.56±4.84f 

0.19±0.97i 


1.35 
-9.85 
1.68 



Table 4: Properties of the non-Gaussian fixed points found in the i? 2 -truncation of 
three-dimensional QEG including all metric fluctuations [24]. 



a valid and systematic approximatation of the underlying fixed function. We stress, 
that these properties are not put in "by hand" . They arise dynamically when solving 
the fixed function equation. 

At this stage, it is illustrative to compare our fixed functions to previous studies 
of the gravitational RG-flow on finite-dimensional subspaces. In [24] the NGFP struc- 
ture arising from the beta-functions of three-dimensional QEG in the i? 2 -truncation 
has been studied in detail. The corresponding flow equations took into account all 
metric fluctuations, but worked with a polynomial expansion up R 2 only. The beta- 
functions on this three-dimensional truncation of the theory space give rise to one 
"physical" NGFP and two "unphysical" NGFPs. Their properties are summarized 
in Table § 

Comparing these results to the ones for conformally reduced QEG reported in 
Tables ^| and [|, we observe systematic differences: the conformal approximation and 
the full computation systematically lead to different signs for the fixed point values 
for Newtons constant and the massive graviton mass m 2 . In addition there is no 
clear relation between the stability coefficients 9 found in the two computations. This 
lends itself to the interpretation that conformally reduced gravity (as implemented in 
this paper) and the full theory may actually probe different conformal field theories. 

In a sense, the results reported in Tables |2| and |3|, are much closer to the properties 
of the matter-induced non-Gaussian fixed points studied in [27,34,67]. Integrating 
out the quantum fluctuations of a minimally coupled scalar field, the NGFP induced 
in the gravitational sector characteristically also appears at negative Newtons con- 
stant g* < 0. This is very much in the spirit of the proposal by 't Hooft [56] that, from 
a path-integral perspective, the quantum fluctuations in the conformal factor should 
be treated as matter fields. Their effect on the running of the gravitational coupling 
constants is qualitatively different from the contribution of the gravitational degrees 
of freedom defining the light-cone-structure of space-time. We hope to further clarify 
this point in the future. 

We close our discussion with a remark on new massive gravity [64,65]. Comparing 
the expansion of our solutions for r = R/k 2 ^> 1 and r < 1, we note that their 
structure is manifestly different. For r <C 1 our fixed functionals indeed have the 
structure of the new massive gravity Lagrangian, while the exact form of the effective 
action is given in (571). We expect that this actually constitutes an interesting 
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observation when discussing unitarity questions in this framework. 

Clearly, it is desirable to generalize our results in various ways. Concerning 
the three-dimensional framework, the next logical step is the inclusion of the tenso- 
rial fluctuations. This would require generalizing the exact asymptotic heat-kernel 
formula ( [2.26 ) to higher order spin fields. This setup would lend itself to a direct 
comparison to the N = 2 polynomial fixed point structure [24] summarized in Table 
Of course, it would also be interesting to revisit the flow of /(i2)-gravity in four 
dimensions [26,49] where, so far, no complete solutions of the fixed function equation 
have been found. Based on the analysis of the RG-flow in the i? 2 -truncation [24], 
we expect that a detailed understanding of the zero-mode contributions originating 
from the transverse-traceless decomposition of the metric fluctuations will constitute 
a crucial prerequisite. We hope to come back to these points in future work. 
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